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1 Introduction 

There are various highly developed methods for establishing semiclassical approxima- 
tions. Probably the most refined method is based on pseudo-differential and Fourier 
integral operator calculi. This extremely technical approach is well suited for getting 
good or even sharp error estimates. Here, sharp refers to the optimal exponent of 
the semiclassical parameter in the error term. These sharp estimates however often 
require strong regularity assumptions on the operators being investigated. 

A different and very simple method based on coherent states gives the leading 
order semiclassical asymptotics under optimal regularity assumptions. 

The method of coherent states was used by Thirring [20] and Lieb [8] to give a very 
short and simple proof of the Thomas-Fermi energy asymptotics of large atoms and 
molecules; see also a recent improvement by Balodis Matesanz and Solovej [16]. This 
asymptotics had been first proved by Lieb and Simon in [7] using a Dirichlet-Neumann 
bracketing method. 

Because of the Coulomb singularity of the atomic potential the pseudo-differential 
techniques are not immediately applicable to the Thomas-Fermi asymptotics. 

In fact, although the Coulomb singularity does not affect the leading order 
Thomas-Fermi asymptotics, in the sense that it is purely semi-classical, it does cause 
the first correction to be of a non-semiclassical nature. The first correction to the 
Thomas-Fermi asymptotics was predicted by Scott in [18] and was later generalized 
to molecules and formulated as a clear mathematical conjecture in [8]. 

The first mathematical proof of the Scott correction for atoms was given by 
Hughes [4] (a lower bound) and by Siedentop and Weikard [19] (both bounds) by 
WKB type methods. Bach [1] proved the Scott correction for ions. 

In [5], Ivrii and Sigal finally managed to apply Fourier integral operator methods to 
the atomic problem and proved the Scott correction for molecules, which was recently 
extended to matter by Balodis Matesanz [15]. 

In [3], Fefferman and Seco gave a rigorous derivation of the next correction (after 
the Scott correction) in the asymptotics of the energy of atoms. This next correction 
had been predicted by Dirac [2] and Schwinger [17]. 

As we shall explain below (see Page 12) one cannot expect to be able to derive 
the Scott correction using the traditional method of coherent states. 

In this paper we introduce a new semiclassical approach generalizing the method 
of coherent states and show that this approach can be used to give a fairly simple 
derivation of the Scott correction for molecules. 

The standard coherent states method is based on representing operators on L 2 (R n ) 
as integrals of the form 




(1) 



where a(u, q) is a function (the symbol of the operator) on the classical phase space 
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M. and H u ,q is a non-negative operator with the properties 

TrU UQ = i, [ n Mq / f M f 9 = 1. 

For the classical coherent states U U!q is the one-dimensional projection \u, q) (u, q\ 
onto the normalized function 

(x\u, q) = (irh)- n/4 e- {x - u)2/2h e iqx/h . (2) 

We generalize this by representing operators in the form 

J GusAusGus ^ nh y ( 3 ) 

Here G u ,q is some self-adjoint operator such that its square plays the role of IT U >q and 
A uq = Bq(u, q) + Bi(u, q) ■ x — ihB 2 (u, q) ■ V is a differential operator linear in x and 
—ihV. (We have denoted by x the position operator.) We shall make an explicit choice 
of Q UA in Sect. 3. In other words, we allow the symbol in the coherent state operator 
representation to be not just a real function on phase space but to take values in first 
order differential operators. If we consider, for example, a Schrodinger operator of the 
form —h 2 A + V(x), where a natural choice of the coherent state symbol would be 
a(u, q) = q 2 + V(u), then the new idea is now to choose the linear approximation 

A u , q = a(u,q) + d u a(u,q)(x - u) + d q a(u,q)(-ihV - q). 

The representation (3) will then be a better approximation of the Schrodinger operator 
than (1) (see Theorem 11 for details). 

In order to explain the Scott correction we consider the non-relativistic Schrodinger 
operator for a neutral molecule 

Z 1 
H(Z, R) = H(Z ± , ...,Z M ;R 1 ,...,R M ) = J2 - V(Z, R, Xi )) + ^ — ■ 

X 7 Ob 4 

i=l i<] 1 ■> 1 

We have \Z\ = ^2^ =1 Zj electrons of charges e = — 1 and M nuclei of charges Zj 
located at the fixed positions Ri, . . . ,Rn', we use atomic units where h 2 = m. The 
number M is an arbitrary but fixed integer throughout this paper. The potential 

M 

V{Z,R,x) = 1 £^ m (4) 

j = l 3 

describes the interaction of a single electron with all the nuclei. The operator H (Z, R) 
acts as an unbounded operator in the space 

A L 2 (R 3 x {-1, 1}), where ±1 refer to 
the spin variables. We are interested in the ground state energy 



E(Z,R) = inispecH(Z,R) 
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and, in particular, in the asymptotic expansion for large charges. In defining the energy 
we have ignored the nuclear repulsion. It would simply shift the energy by a constant 
depending on Z and R. 

The version of the Scott correction that we prove in this paper can now be stated 
as follows. 

Theorem 1 (Scott correction). Let Z = \Z\(z 1: . . . , z M ), where zi, . . . , z M > and 

R = |Z| _1//3 (ri, . . . ,tm), where |rj — rj\ > r for some r > 0. Then, 

E(Z,R) = E TF (Z,R) + l Z] + 0(\Z\ 2 - 1 / 30 ), (5) 

l<j<M 

as \Z\ — > oo ; where the error term (9(|Z| 2-1 / 30 ) besides \Z\ depends only ore ^i, ... , zm, 
and r . 

This is established in lemmas 18 and 19. In fact, one could improve slightly on the 
error estimate to the expense of limiting the range of Z and R, and vice verse. 

It turns out that E TF (Z,R) is of order \Z\ 7 ^ 3 and the next term | Xa<j<M -^j * s 
the Scott correction. 

Part of our derivation of Theorem 1 is similar to the multi-scale analysis in [5] 
and we adopt their notation. Our semiclassical method, however, is very different. It 
does not rely on the spectral calculus, but uses only the quadratic form representa- 
tion of operators. Moreover, we treat the Coulomb singularities completely differently 
from [5]. In treating the singularities and the region near infinity the Lieb-Thirring 
inequality plays an essential role. 

Another virtue of our proof is that it gives an explicit trial state for the energy 
that is correct to an order including the Scott correction. This is, in fact, how we 
prove that the Scott correction is correct as an asymptotic upper bound. 

This paper is organized as follows. In Sect. 2.1 we list for the convenience of the 
reader the analytic tools that we shall use in a crucial way. In Sect. 2.2 we review 
Thomas- Fermi theory. In Sect. 3 we introduce the new coherent states. In Sect. 4 we 
apply this new tool to prove the semi-classical expansion of the sum of the negative 
eigenvalues of a non-singular Schrodinger operator localized in some bounded region of 
space. This is the key application of our new method. The proof for the semi-classical 
expansion for the Thomas- Fermi potential is presented in Sect. 5. In Sect. 6 we finally 
prove lower and upper bound for the molecular quantum ground state energy. Some 
calculations concerning the new coherent states and a theorem on constructing a 
particular partition of unity are put into the appendices. 

2 Preliminaries 

2.1 Analytic tools 

In this subsection we collect the main analytic tools which we shall use throughout the 
paper. We do not prove them here but give the standard references. Various constants 
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are typically denoted by the same letter C, although their value might, for instance, 
change from one to the next line. 

Let p > 1, then a complex-valued function / (and only those will be considered 
here) is said to be in L p (W n ) if the norm ||/|| p := (J \f(x)\ p dx) 1 ^ is finite. For any 
l<p<t<q<oowe have the inclusion L p H L q C L*, since by Holder's inequality 

11/11* < ll/llp A ll/lir A with X?- 1 + (1 - A)^- 1 = t~\ 

We call 7 a density matrix on L 2 (R") if it is a trace class operator on L 2 (IR n ) 
satisfying the operator inequality < 7 < 1. The density of a density matrix 7 is the 
L 1 function p 7 such that Tr(7#) = J p 1 (x)9(x)dx for all 9 e C£°(M. n ) considered as 
multiplication operators. 

We also need an extension to many-particle states. Let ip G (§Z) N L 2 (R 3 x { — 1,1}) 
be an iV-body wave-function. Its one-particle density p^ is defined by 

N 

PA X ) = ^2 ^ • • • ^ / l^/K^i? s i • • • , x n, s N )\ 2 5(xi - x) dxi ■ ■ -x N . 

i=l si=±l Sjv=±1 "* 

The next inequality we recall is crucial to most of our estimates. 

Theorem 2 (Lieb-Thirring inequality). One-body case: Let 7 be a density op- 
erator on L 2 (W l ), then we have the Lieb-Thirring inequality 

Tr[-|A 7 ] >K n j p\^ n (6) 

with some positive constant K n . Equivalently, let V E L 1+ "/ 2 (lR n ) and 7 a density 
operator, then 

Tr[(-|A + y) 7 ] > -L n J |\/_| 1+ "/ 2 , (7) 

where x_ := min{x, 0} 7 and L n some positive constant. 

Many-body case: Let ip e /\ N L 2 (R 3 x {-1, 1}). Then, 

(Vf>±A^> 2-*<*K 3 Jpf. (8) 

The original proofs of these inequalities can be found in [6]. 

From the min-max principle it is clear that the right side of (7) is in fact a lower 
bound on the sum of the negative eigenvalues of the operator — |A + V . 
We shall use the following standard notation: 

£>(/) = £>(/,/) = ! J J f(x)\x-y\- 1 f(y)dxdy. 
It is not difficult to see (by Fourier transformation) that ||/|| := D(f) 1 ^ 2 is a norm. 
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Theorem 3 (Hardy-Littlewood-Sobolev inequality). There exists a constant C 
such that 

D(f)<C\\f\\l /5 . (9) 

The sharp constant C has been found by Lieb [11], see also [12]. 
In order to localize into different regions of space we shall use the standard IMS- 
formula 

-\9 2 A - \M 2 = -9A9 - (V0) 2 , (10) 

which holds, by a straightforward calculation, for all bounded C 1 -functions 9 (here 
considered as a multiplication operator). 

Finally we state the two inequalities which we need to estimate the many-body 
ground state energy E(Z, R) by an energy of an effective one-particle quantum system. 
The first one is an electrostatic inequality providing us with a lower bound. This 
inequality is due to Lieb [10], and was improved in [13]. 

Theorem 4 (Lieb-Oxford inequality). Let ip E L 2 (R m ) be normalized, and p^ its 
one-electron density. Then, 




> D{ P ,,) - C j ' pf. (11) 



The best-known numerical value is 1.68, but this does not play a role here. 
An upper bound to E(Z, R) is provided by a variational principle for Fermionic 
systems. This is also due to Lieb [9] . 

Theorem 5 (Lieb's Variational Principle). Let 7 be a density matrix on L 2 (R 3 ) 
satisfying 2Tvy = 2 f p^{x) dx < Z (i.e., less than or equal to the number of electrons) 
with the kernel p-y(x) = ^(x,x). Then 

E{Z,R) <2Tr\(-^A-V(Z,R,x))j\ + D{2 Pl ). (12) 
The factors 2 above are due to the spin degeneracy. 



2.2 Thomas- Fermi theory 

Consider z = (z u . . . , z M ) G R+ and r = (n, . . . , r M ) G M 3M . Let < p G L 5 / 3 (R 3 ) n 
L^IR 3 ) then the Thomas- Fermi (TF) energy functional, £ TF , is defined as 

£ TF (P) = ^(3vr 2 ) 2/3 J p{xfl* dx- J V(z, r, x)p(x) dx + D(p), 

where V is as in (4). 

By the Hardy-Littlewood-Sobolev inequality the Coulomb energy, D(p), is finite 
for functions p E L 5 / 3 (M 3 ) n L 1 (R 3 ) C L 6 / 5 (R 3 ). Therefore, the TF-energy functional 
is well-defined. Here we only state the properties about TF-theory which we use 
throughout the paper without proving them. The original proofs can be found in [7] 
and [8]. 
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Theorem 6 (Thomas-Fermi minimizer). For all z = (zi,...,Zm) £ o.nd 
r = (ri, . . . ,tm) £ K 3M toere exzsfo a unique non-negative p TF (z,r,x) such that 
f p TF (z, r, x) dx = J2kLi z k an d 

^ TF(p TF ) = inf |^TF (p) . o^ei^tfjni 1 ^)}. 

We s/ia// denote by E TF (z,r) := £ TF (p TF ) i/ie TF-energy. Moreover, let 

V TF (z, r, x) := V(z, r, x) - p TF * |a;| _1 . (13) 

fre the TF -potential, then V TF > and p TF > ; and p TF is i/ie unique solution in 
L 5 / 3 (R 3 ) n L 1 (R 3 ) to toe TF-equation: 

V TF (z, r, x) = \ (3tt 2 ) 2 V F (z, r, x) 2 / 3 . (14) 

Very crucial for a semi-classical approach is the scaling behavior of the TF- 
potential. It says that for any positive parameter h 



V TF (z,r,x) = /i-V TF (/i 3 z,/i- 1 r,/i- 1 a;), (15) 
p TF (z,r,a;) = h~ 6 p TF '(/i 3 z, /rt, /TV) (16) 
£ TF (z,r) = /i- 7 £ TF (/i 3 z,/rt). (17) 

By h~ l r we mean that each coordinate is scaled by h' 1 , and likewise for h 3 z. By the 
TF-equation (14), the equations (15) and (16) are obviously equivalent. Notice that 
the Coulomb-potential, V, has the claimed scaling behavior. The rest follows from the 
uniqueness of the solution of the TF-energy functional. 

We shall now establish the crucial estimates that we need about the TF potential. 

Let 

d(x) = mm{\x - r k \ \ k = 1, . . . , M} (18) 

and 

f(x) = min^x)- 1 / 2 , d(x)' 2 }. (19) 
For each k— 1, . . . , M we define the function 

W k (z,r,x) = V TF (z,r,x)-z k \x-r k \- 1 . (20) 

The function W k can be continuously extended to x = r k . 

The first estimate in the next theorem is very similar to a corresponding estimate 
in [5]. 

Theorem 7 (Estimate on V TF ). Let z = (z 1: . . . , z M ) £ and r = (r 1: . . . , r M ) £ 
R 3M . For all multi-indices a and all x with d(x) ^ we have 

|^\/ TF (z,r,x)| < C a f(x) 2 d(x)-W, (21) 
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where C a > is a constant which depends on a, z±, . . . ,zm, and M. 

Moreover, for \x — r k \ < r min /2 7 where r min = min^ \r k — re\ we have 

0<W k (z,r,x)<Cr-^ + C, (22) 

where the constants C > here depend on zi, . . . , zm , and M . 

Proof. Throughout the proof we shall denote all constants that depend on a, 
zi, . . . , zm, M by C a . Constants that depend on zi, . . . , zm we denote by C. In this 
proof we shall omit the dependence on r and z and simply write V TF (x) and Wk(x). 
We proceed by induction over |a|. If a — we have the well known bound [7] that 

M 

< max{Kf (x) | k= 1, . . . , M} < V TF (x) < £ V™(x), (23) 

fe=i 

where V^ F denotes the Thomas-Fermi potential of a neutral atom with a nucleus 
placed at r k e ffi 3 with nuclear charge z k . This potential satisfies the bounds [7] 

C_ min{z fc |x - r fc | _1 |x - r fc |" 4 } < \/^ F < C+ min{^|a; - r fe | _1 , |x - r fe |" 4 }, (24) 

where C± > are universal constants (note that by scaling (15) it is enough to 
consider the case z^ — l). We therefore get that 

C_ min{^, . . . , z M , l}f(x) 2 < V TF (x) < C + M max{^, ...,z M , l}f(x) 2 . (25) 

This in particular gives (21) for a — 0. Assume now that (21) has been proved for all 
multi-indices a with \a\ < M, for some M > 0. We shall first establish an estimate 
for the derivatives d a p of the TF density p. 

From the TF equation we have that p = C(V TF ) 3 ^ 2 . Thus d a p(x) is a sum of terms 
of the form 

V TF (x) 3 / 2 - k d^V TF (x) ■ ■■d^ F (x) 

where k — 0, . . . , |a| and \f3\\ + . . . + \f3k\ = \a\. Thus by the induction hypothesis and 
(25) we have for |a| < M that 

\d a p(x)\ <C a f{xfd{x)-\ a y (26) 

We now turn to the potential. Given a with \a\ = M. Choose some decomposition 
a — (5 + a', where \/3\ = 1 and \a'\ = M — 1. 
For all y such that \y — x\ < d{x)/2 we write 

d a 'V TF (y) = - [ d a 'p{u)\y - u\- l du + R(y), 

J\u-x\<d(x)/2 

where 

R(y) = d a 'V TF (y) + / d a 'p(u)\y - u\~ l du 

J\u-x\<d(x)/2 
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is a harmonic function of y for \y — x\ < d(x)/2, since AV TF (y) = A.irp{y) for such y. 
It follows that 

\d p R{x)\ < Cdix)' 1 sup 



|£-s|=d(x)/2 



This can be seen by the Poisson formula 

d(xf 



R(y) = j 



-\y-x\ 



\£-x\=d{x)/2 



27rd(x)\Z-y\ 3 



valid for \y — x\ < d(x)/2. Here dS denotes the surface measure of {£ : |£ — x\ — 
d(x)/2}. To estimate sup,^_ x < d ^ x y 2 \R(£)\ we use the induction hypothesis, i.e., (21) 
and (26) for a'. Note that for |f - x\ < d(x)/2 we have that d(x)/2 < d(t) < 3d(x)/2 
and hence also that /({) < 4/(x). Thus, if we also use that f(x) < d(x)~ 2 we get that 

sup \R{i)\<C a ,d{x)-\ a '\f{x) 2 

\a- x \=d{x)/2 



and hence that 



\dPR(x)\ <C a d{x)-\ a \f{xf. 



Finally we have that 

a?/ 

J\u 

The estimate (21) follows since 



d a ' p{u) \y — u\ l du 


<-j 


d a 'p(u) 


\y-u\ 


x\<d(x)/2 


J\u-x\<d(x)/2 







I 



\u-x\<d(x)/2 



d a ' p{u) 



\X — U\ 



2 du < C a ,f(xfd(x)-\ a '\ +l < C a f(xfd(x)-\ a \ 



where we have again used that f(x) < d(x)~ 2 . 

The estimate (22) follows from (23) and (24) if we note that the atomic potential 
satisfies 

< V r Y(x) - z k \x - r k \~ l = p T * * \x\~ l < C. 
Here p^f is the atomic density. These last estimates follow since pJ F is non-negative 



and bounded in L 5 / 3 and in L 1 . 



□ 



The relation of Thomas-Fermi theory to semiclassical analysis is that the semi- 
classical density of a gas of non-interacting electrons moving in the Thomas-Fermi 
potential V TF is simply the Thomas-Fermi density. More precisely, the semiclassical 
approximation to the density of the projection onto the eigenspace corresponding to 
negative eigenvalues of the Hamiltonian — |A — V TF is 



ly 



l 



. 2 -V TF (z,r,:r)<0 



dp 
(2^ 



2 3 / 2 (3tt 2 )- 1 (\/ tf ) 3 / 2 (z, r, x) = p TF (z, r, x). 
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Here the factor two on the very left is due to the spin degeneracy. Similarly, the 
semiclassical approximation to the energy of the gas, i.e., to the sum of the negative 
eigenvalues of — |A — V TF is 

2 / (i? 2 - V™(z, r, *)) ^ = -^/ V' TF (z, r, = £™(z, r)+B(„ TF ). 

(27) 

In Section 5 we shall make the semiclassical approximation more precise. 



3 The new coherent states 

We shall now define the new coherent states (or better coherent operators 1 ) discussed 
in the introduction. I.e., we shall define the operator Q uq that is used to represent 
operators in the form (3). 

The classical coherent states (2) localize in both x and p on a scale of order h. We 
shall choose Q Ujq to localize on a longer scale. We define 



(28) 



The new scale is 1/a > h. Note that if we let a — > 1/h then Q 2 q converges to the 
projection n Ui9 = \u,q)(u,q\. A straightforward calculation gives the following result. 

Lemma 8 (Completeness of new coherent states). These new coherent operators 
satisfy 

I (2^ = G '<* " «>• / 6 « wk» = Gt( - !ftV " ql 

where x denotes the operator multiplication by the position variable x. Here G^v) = 
(b/ir) n / 2 exp(— bv 2 ) with b = 2a/(l + h 2 a 2 ). Note that Gb has integral 1 and hence 



I 



2 dudq 
q {2nh) n 



We shall study operators that can be written in the form (3). If A U:q = B (u, q) + 
B\(u, q) •x — ihB2{u, q) -V is the operator valued symbol in (3) we shall denote by A UA 
the linear function A UA {y,p) = B (u, q) + B^u, q) • v + B 2 (u, q) ■ p. When A u q (v,p) 
is independent of (v,p), i.e., if B\ = B 2 = and if a — > h^ 1 we recover the usual 
coherent states representation of an operator. Thus, on the one hand we do not use 
as sharp a phase space localization as the one-dimensional coherent state projection 
since a < 1/h, but on the other hand, we use a better approximation than if A u ^ q were 
just a constant. 

1 This should not be confused with the quantum coherent operators introduced by Lieb and Solovej 
in [14] in order to compare two quantum systems. 
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More generally we shall consider operators of the form 

Qu, q f (A u , q ) G u ,qdudq, (29) 



/' 



where / : R — > R is any polynomially bounded real function. As we shall see in 
the next theorem the integrand above is a traceclass operator for each (u,q). The 
integral above is to be understood in the weak sense, i.e., as a quadratic form. We 
shall consider situations where the integral defines bounded or unbounded operators. 

Theorem 9 (Trace identity). Let f : R -> R and V : R n -> R be polynomially 
bounded, real measurable functions and 

A = B + B x x - ihB 2 V 

a first order self-adjoint differential operator 2 with B E R, B 12 E R n . Then 
Q u>q f(A) Q u q V (x) is a trace class operator (when extended from C£°(M. n )) and 

Tr[g u , q f(A) Q UA V(x)} = J f(B + B lV + B 2 p) G b {v - u)G b (q - p) 

x G(p h 2yi (z)V(v + h 2 ab{u — v) + z)dvdpdz. 



In particular, Tr 



The proof is given in Appendix A. We shall need the following extension of this 
theorem, where we however only give an estimate on the trace. The proof is again 
deferred to Appendix A. 

Theorem 10 (Trace estimates). Let f,A be as in the previous theorem. Let more- 
over <p E C n+4 (R n ) be a bounded, real function with all derivatives up to order n + 4 
bounded and V,Fe C 2 (M. n ) be real functions with bounded second derivatives. Then, 
for h < 1, 1 < a < 1/h and b = 2a/(l + h 2 a 2 ) we have with a(u,q) = F(q) + V{u) 
that 3 



'71 o 

/ 



X 



Tr[g Utq f(A)g u , q <j>(x) (F(-ihV)+V(x))<f>(x)\ 
f{B + B lV + B 2 p) G b {v - u)G b {q - p) 

(v + h 2 ab{u — v)) 2 + Ei(u, v]^o-(v + h 2 ab{u — v),p + h 2 ab{q — p)) 

+ E 2 (u,v;q,p) dvdp, 

with \\Ei\loo, 1 1 £2 1 |oo < Ch 2 b where C depends only on 

sup Hd^Hoo, sup Hd'VHoo, and sup H^FHoq. 

\u\<n+4 \v\=2 |i/|=2 

(Note that the assumption 1 < a < 1/h implies 1 < b < 1/h.) 
2 The operator A is essentially self-adjoint on Schwartz functions on W 1 . 

3 The operator G u , q f(A) G u , g (f>(x) (F(-ihV) + V{x))<t>{x) is originally defined on, say C °°(R 11 ), 
but it is part of the claim of the theorem that it extends to a traceclass operator on all of L 2 (R"). 
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The above theorem shall be used to prove an upper bound on the sum of eigen- 
values of the operator F(—ihV) + V(x), in the case when F(p) = p 2 . This is done in 
Lemma 14 by constructing a trial density matrix in the form (29). 

To prove a lower bound on the sum of the negative eigenvalues one approximates 
the Hamiltonian F(—iKS7) + V(x) by an operator represented in the form (3). This 
approximation which we now formulate is also proved in Appendix A. 

Theorem 11 (Coherent states representation). Consider functions F,V E 
C 3 (R n ) ; for which all second and third derivatives are bounded. Let o~(u,q) = F(q) + 
V(u), then we have for a < 1/h and b = 2a/(l + h 2 a 2 ) the representation 

F(-ihW) + V(x) = J Q u , q H u , q Gu, q -^ + E 

(as quadratic forms on C£°(M. n )), with the operator-valued symbol 

H u , q = a(u, q) + ^-Aa(u, q) + d u a(u, q)(x - u) + d q a(u, q)(-ihV - q). (30) 
The error term, E is a bounded operator with 

l|E|| < cr 3/2 Halloo + ch 2 b \\d a o-\\oo- 

H=3 \a\=2 



4 Semi-classical estimate of Tr <f>(—h 2 A + V)(f\ 

In this Section we study the sum of the negative eigenvalues of Schrodinger operators 
with regular potentials localized in a bounded region by a localization function 0. 
This will turn out to be the key theorem in this paper. Let us recall our convention 
that X- — (x)- — min{a;, 0}. For convenience we consider balls, and we start with the 
unit ball. 

Theorem 12 (Local semiclassics). Let n > 3 and <fi e CQ +4 (lR n ) 7 be supported in 
a ball B dW 1 of radius 1 and V G C 3 (B) be a real function. Let H = -h 2 A + V, 
h > and a(u, q) = q 2 + V(u) . Then, 



Tr[0#0]_ - (27rh)~ n J <p 2 (u)o-(u,q)^dudq 



< Ch- n+%/ \ 



The constant C > here depends only on n, \\(fi\\c n + 4 and 1 1 1 1 c 3 - [Here ||T^||c 3 — 

BUP|a|<sll^lloo.y 

With the classical coherent states the estimate one would normally prove would 
be that the right side above is Ch~ n+1 . We find it instructive to sketch the proof of 
this here in order to make the comparison with the new method clearer. 
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We may assume that V is defined on all of R n with bounded second and third 
order derivatives. We shall here assume that h < 1. From Theorem 11 with a — 1/h 
we have 

dudq 

+ E, 



H = -h 2 A + V = J[q 2 + V(u)} \u, q) (u, q\ 



{2iih) n 

with ||E|| < Ch. The constant depends on the second and third order derivatives of 
V, which are bounded. We have here used that for a = b = 1/h the first order terms 
in x and V do not contribute (see (99) below). Moreover, the term Aa is of order h 
since V has bounded second order derivatives. The error term E can be controlled 
using the Lieb-Thirring inequality as in (31) below. 
Then from Theorem 10 we have 



Tr[0if0]_ > Jiq 2 + V(u)]_Tr[<P\u,q)(u,q\ ( p 



/ 



(2nh) n 
dudq ^ L _ n+ i 



dudq -Ch- n+1 



[ q 2 + V(u)]^\u) WJ ^-Ch 



where C now also depends on the derivatives of 0. 
For an upper bound we set 



l=j X(-oo,o][g 2 + V(u)} \u,q)(u,q\ 



(2tt/i)"' 



where X(-oo,o] denotes the characteristic function of the interval (— oo, 0]. It is clear 
that < 7 < 1. When calculating Tr[70if0] we may again refer to the general 
theorem 10 with a = b = 1/h. We obtain 



X(-oo,o][<? 2 + V( u )] Tr \u, q){u, q\4>H(f) 



. (2vr/i) n 



As mentioned in the Introduction it is important that we obtain errors bounded 
by Ch~ n+1+£ for some e > as in Theorem 12. We shall prove Theorem 12 by again 
proving upper and lower bounds on Tr[0iJ0]_. 

Lemma 13 (Lower bound on Tr(0if0)_). Let n > 3, <fi E Cq +4 (W 1 ) be supported 
in a ball B of radius 1 and assume that V G C 3 (B). Let H = —h 2 A + V, h > 0. Then, 

Tr[0#0]_ > (2nh)- n J <p 2 (u)a(u, q)^dudq - Chr n ^ jh . 



The constant C > here depends only on n, ||0||o+ 4 an d \\V\ 



Proof. Since has support in the ball B we may without loss of generality assume 
that V G Cq(R 3 ) with the support in a ball B 2 of radius 2 and that the norm ||T^||c 3 
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refers to the supremum over all of M. n . We shall not explicitly follow how the error 
terms depend on ||0||c 3 an d ||^||c 3 - All constants denoted by C depend on n, ||0||c 3 > 



c 3 - 



First note that by the Lieb-Thirring inequality we have that 
Tr[0#0]_ > CHWl 



*(u, q )_*%->-Ch-«. 

Consider some fixed < r < 1 (independent of h). If h > r then 

Tr[0#0]_> J ^u)a(u, q )_^^-Cr-^h-^. 

We are therefore left with considering h < r. Of course one should really try to find 
the optimal value of r (depending on 0, and V) we shall however not do that. In 
studying the case h < r it will be necessary to assume that the choice of r is small 
enough. We therefore now assume that h < r and that r is small. 
From Theorem 11 we have that 



Tr[0#0]_ > Tr 



1 G u,qH u,q@ u,q(i 



dudq 



(2nhy 

+Tr [0 (-eh 2 A - C{b^ 2 + h 2 b)) 0] 



(31) 



where < e < 1/2 and 



H u , q = <t(u, q) + — Aa(u, q) + d u a(u, q)(x - u) + d q a(u, q)(-ihV - q) 

with a(u, q) = (1 — e)q 2 + V(u). We shall choose a depending on h satisfying r _1 < 
a < h^ 1 and hence r _1 < b < h^ 1 . It is clear (e.g. from the Lieb-Thirring inequality) 
that the second trace above is estimated below by —Ch~ n e~ n / 2 (b~ 3 / 2 + h 2 b) l+n / 2 . We 
shall choose e = \(b~ 3 ^ 2 + h 2 b); note that e < 1/2. Thus we find that the second trace 
is estimated by —Ch^ n (b^/ 2 + h 2 b). From the variational principle we have 



Tr[0#0]_ > / Tr 



4>Gu,q 



H, 



"A 



Gu,q4> 



dudq 



-Ch- n (b- 3/2 + h 2 b). 



We first consider the integral over u outside the ball B 2 of radius 2, where V = 0. 
Using Theorem 9 (with V replaced by 2 ) and J 4> 2 < C, we get that this part of the 
integral is 



/ 



n 



1 " £ )<T + 7n(l - e) + 2(1 - e)q ■ (p - q) G b (p - q)G b (u - v) 



2b 



x G(ph2}-i(z)(f)(v + h 2 ab(u — v) + z) 2 dvdpdz 



dudq 
(2nh) n 



> CJ(1-S)[ P 2 -( P - q) 2 ]^G b {p - q)^L > _ C76 -(n + 2)/2 /l -„ j 
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which for all dimensions n is bounded below by —Cb~ z l 2 hr n . Actually it is not difficult 
to see that we could have inserted a factor e~ Cb on the right of this estimate since 
u B 2 and is supported in Bi, but we do not need this here. 

For the integral over u G B 2 we use Theorem 10 with F = and V — 1 to obtain 

Tr[0#0]_ > I U (v + h 2 ab(u - v)) 2 + Ch 2 b) G b {u - v)G b {q - p) 

x [H u , q (v,p)}_ ^^dpdv - Ch~ n (b-V 2 + h 2 b), (32) 

where 

H u , q (v,p) = °(u, q) + ^ A3r ( M ' + d u a(u, q)(v - u) + Oqcr(u, q)(p - g). 

The rest of the proof is simply an estimate of the integral in (32). Note that by Taylor's 
formula for a we have 

H u , q (v,p) > v(v,p) +£ v (u-v,q-p) -C\u-v\{b~ l + \u-v\ 2 ), (33) 

where 

q) = ^ A5r ( w ' °) - ( x - £ )<7 2 ~ 2 S ^^(v)^-. 

We have here used that Aa(v,p) is independent of p and that \Aa(v, 0) — Aa(u, 0)| < 
C|it — v\. Since ||U||c3 < oo and thus, in particular, a{v ,p) > (1 — e)p 2 — C we easily 
get that 

' G b {u - v )G b {q-p) [H u , q (v,p)}_ dpdqdv > -C (34) 



and hence from (32) that 

Tr[0#0]_ > I <f)(v + h 2 ab(u-v)) 2 G b (u-v)G b (q-p) 

x [H u , q (v,p)]_ ^^dpdv - Ch~ n (b-'/ 2 + h 2 b). 

Here we have of course used the fact that the ^-integration is over a bounded region. 
From now on we may however ignore the restriction on the -u-integration. Using (33) 
we find after the simple change of variables u — > u + v and q — > q + p that 

Tr[0#0]_ > J 4> (v + h 2 abu) 2 G b {u)G b {q) 

dudq 



a(v,p) + £ v (u,q) - C\u\(b 1 + \u\ 
Ch- n (b-*l 2 + h 2 b). 



(2nh)' 



dpdv 
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We now perform the p-integration explicitly. Recall that a(v,p) = (1 — e)p 2 + V(v) 
and that J (p 2 + s)~dp = — ^|^c<j n |s_|(™/ 2 ) +1 , where oo n is the volume of the unit ball 
in lR n . We get 

Tr[0#0]_ > — ^(l- e )-5 [ <f) (v + h 2 abu) 2 G b {u)G b {q) 
n + 2 J 



x 



V{v) + Uu,q) ~ CluKb- 1 + \u\ 2 ) 



'r + l 



dudq 



dv 



-Ch- n (b-^ 2 + h 2 b). 
By expanding we find that 



(2nh) n 
(35) 



-+1 

2 ^ 



< 



\V(v)^ +1 - (^ + i)\v(vU n H v (u, q ) 



+c 



Uu,q) +C\u\(b- 1 + \u\ 2 ) + \u\ 2 ). 

We have here used that since n > 3, the function K9ih is C 2 . Hence 

2uj r , 



Tr[0F0]_ > -_^L(l- £ )-f J ct>{v + h 2 abu) 2 G b {u)G h {q) 



dudq 
(2nh) n 



dv 



< Ch 4 a 2 b 2 \u\ 2 < CW\u\\ 



2i.2l.i2 



x (jv>)_|^-(_ + l)|V>)^Kg) 

-C7r n (r 3 / 2 + /^). 

We now expand 4> 2 

<f)(v + h 2 abu) 2 - <f){v) 2 - h 2 abu ■ V (0 2 ) (u) 
and use the crucial identities 

J £ v (u, q)Gb{u)Gb{q)dudq = and ^ uGb(u)du = 0. 

We thus arrive at 

(27r/i)"Tr[0#0]_ > -^-(l - J (t>(v) 2 \V(v)^ +1 dv - C(b~^ 2 + h 2 b) 
= (l-e)"5 J (f ) (v) 2 a(v 1 p)_dvdp-C(b- z/2 + h 2 b). 

The lemma follows if we choose a = max{/i _4//5 , r -1 } and £ = i(6~ 3 / 2 + /i 2 6). Recall 
that a < b < 2a. Thus b^' 2 < a^l 2 < h 6 / 5 and h 2 b < 2h 2 a < 2r~ 1 / 5 h 6 / 5 . 

□ 
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In order to prove an upper bound on Tr(0ff0)_ we shall use that for any density 
matrix 7 (i.e., a traceclass operator with < 7 < 1) we have from the variational 
principle that Tr(0if0)_ < Tr ((j)H(jrf). Hence the upper bound needed to prove The- 
orem 12 is a consequence of the following lemma. 

Lemma 14 (Construction of trial density matrix). Let n > 3, G CQ +4 (IR n ) 
be supported in a ball B of radius 1, and V G C 3 (B). Let H = —h 2 A + V , h > and 
a(u, q) = q 2 + V(u). Then there exists a density matrix 7 on L 2 (M n ) such that 

Tr[0#0 7 ] < J cf> 2 (u)a(u, q)-j^ + Ch~^\ (36) 
Moreover, the density of '7 satisfies 

Pl {x) - (2irh)- n cu n \V(x)_\ n/2 < C/i- n+9/10 , (37) 
for (almost) all x G B and 



J 4>{x) 2 Pl {x)dx - {2nh)- n u) n J <f)(x) 2 \V(x)_\ n/2 dx 



< Ch~ n+6 / 5 , (38) 



where u n is the volume of the unit ball in W 1 . The constants C > in the above 
estimates depend only on n, ||0||c™+ 4 ; an d ll^llc 3 - 



Proof. As in the lower bound we choose some fixed < r < 1. We have for h > t 
that for some C > 



/ 



^u)a(u,q).^^ + Cr-^h-^>0 



and 

\(2nh)- n u; n V(x)_\ n/2 < Cr- 6 / 5 h~ n+6/5 , 

If h > t we may therefore choose 7 = 0. We may therefore now assume that h < r 
and if necessary that r is small enough depending only on 0, and V. Also as in the 
lower bound we may assume that V G Co(lR n ) with support in the ball B 3 / 2 concentric 
with B and of radius 3/2. 

In analogy to the previous proof for the lower bound we now for each (u, q) define 
an operator h u q by 



hu,q 



a(u, q) + ±Aa(u, q) + d u a(u, q)(x - u) + d q a(u, q)(-ihV - q) , if u G B 2 

, if u £ B 2 ' 



The corresponding function is 

j a(u,q) + j E Aa(u,q) + d u a(u,q)(v-u) + d q a(u,q)(p-q) , if u G 5 2 
'VA^W-i ,if^^5 2 
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Recall that b = 2a/ (1 + h 2 a 2 ) (i.e., in particular a < b < 2a) and as in the lower 
bound we shall choose a = max{/i~ 4 / 5 , r -1 } 
Similar to (33) we have for u G B 2 that 



where 



\h u , q (v,p)-a(v,p)-£; v (u-v,q-p)\<C\u-v\(b ^lu-v] 2 ), (39) 
Uu, q) = ^A<7(u, 0)-q 2 -lYl didjVW 



WiUj. 



We have here used that Aa(v,p) is independent of p. 

If we let x — X(-oo,o] be the characteristic function of (— oo, 0] we now define 



7 



Qu,qX[K,q] Q 



u,q 



dudq 
(2nh) n ' 



(40) 



Since < x[hu,q] < 1 it is obvious that < 7 < 1. Moreover, by Theorem 9 and (39), 
7 is easily seen to be a traceclass operator with density 



X {h u , q {v, p)) G b (u - v)G b (p - q)G {bh 2 ) -i (x-v - h 2 ab(u - v))dvdp-^ 



dudq 



h) n ' 

If we change variables u ^ u + v, g — > q + p and perform the p-integration we find 
that 



Pl (x) 



UJr, 



dudq 

E(v, u, q)G b (u)G b (q)G {bh 2)-i (x - v - h 2 abu)dv 



E(v — h 2 abu, u, q)G b (u)G b (q)G^h 2 )- 1 ( x — v)dv 



dudq 
(2nh) n 



(41) 



(l-h 2 ab)ueB 2 -v 

where E(v , u, q) = u^ 1 f x(h(u+v,q+p)(v , p)) dp > 0. From equation (39) we have 

E(v, u, q) 2 > n - J (V(v) + & (u, g)) _ 1 1 < C\u\ (b' 1 + \u\ 2 ), 

for all v , q G 1R™ and u G B 2 — v. Since 

\£ v (u,q) -t,v-h*abu{u,q)\ < Ch 2 ab\u\(b~ l + \u\ 2 ) 
we therefore also have 



(42) 



E(v - h 2 abu, u, q) 2/n - (v(v) + r] v (u, g)) 



< Ch A a 2 b 2 \u\ 2 

+ C(l + h 2 ab)\u\(b- 1 + \u\ 2 ), 



where 



rj v (u, q) = £ v (u, q) — h 2 abW(v) ■ u. 
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Hence from (41) 



V(v) +Vv(u,q)j G b (u)G b (q)G {bh2) -i(x - v)dv 



dudq 



(l-h 2 ab)ueB 2 -v 



< Ch- 2 (h 4 a 2 b + b' 3 / 2 ) < C/ i - 2+6 / 5 ,(43) 



where C may depend on r. 

We now use that for all x, y e M and all n > 3 we have 

n i= i i — i ni i = N I r C\x — n = 3 /..s 

b_ 2 - |/_ 2 + S y_ 2 l (x-y) \ < < ,n ' . ,n' ' |9 ^ , (44) 

" 1 |y 1 2iy 1 v y>] - \ C(\x\2~ 2 + \y\2- 2 )\ x -y\ 2 , n>A y ' 

where C depends on n. This gives for n = 3 (it is left to the reader to write down the 
estimates for n > 4) 



3 

V(v) +rj v (u : q)^_ 2 - \V(v)_\i + §|V»_|^( M ,g) 



< C\r) v (u,q)\*. (45) 



It is now again crucial that J rj v (u,q)G b (u)G b (q)dudq = and hence for v e 

SUpp(V) C .03/2 

y r] v (u,q)G b {u)G b (q)dudq < Ce^ 5 < C/i 6/5 . (46) 

(l~h 2 ab)u£B 2 ~v 

Combining (43), (45), (46), and \r) v (u,q)\ < C(b~ 1 + \u\ 2 + \q\ 2 + h 2 ab\u\) we obtain 

Pl {x) - (2nh)- 3 L; 3 J \V(v)^ 2 G {bh 2 ) - 1 (x-v)dv 

< Ch-\e- h ' 5 + h 3 a 3 ' 2 b 3 / A + b- 3 ' 2 + h 6 ' 5 ) < Ch- 3+& /\ (47) 

where we have again removed the condition (1 — h 2 ab)u G B 2 — v paying a price of 
Ch~ 3 e- h l\ 

A simple Taylor expansion of <fi 2 gives 



4>(x) 



J <p(v) 2 G( bh 2)-i(x - v)dv 



< Cbh 2 < 



where we have again used that J vG( b h2)-i(v)dv = 0. This immediately gives (38). 
Finally, using again (44) we get 

\V(x + v)^ - \V(x)-\* + l\V(x)-\*VV(x)-v < C{\v\i + \v\ 2 ), 

and hence from (47) 

Pl (x) - (2nhT 3 uj n \V(x)^ 2 < Ch' 3 (h 6 / 5 + (bh 2 f 4 ) < Ch- 3+9 / w . 
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We must now calculate Tr(70iJ0) = Tr(70(— h 2 A)(f)) + Tr( / y<j)V<j)) for n > 3. From 
the argument leading to (38) we have 

(27r/i) n Tr( 7 0\/0) < -u n J <j)(x) 2 \V(x)^ +1 dx + Ch- n+6 ^. (48) 

From Theorem 10 we have 

(27T/ i ) n Tr( 7 0(-/ i 2 A)0) = J x(h u , q (v,P))G b (u-v)G b (q-p)[E 2 + 

(<j>(v + h 2 ab(u - v)f + E{){p + h 2 ab(q - p)) 2 ]dudqdvdp, 

where E 1 ,E 2 are functions such that ||£'i|| 0O , H-Eslloo < C/i 2 6. Since 

x(h u ,q(v,p)) G b (u - v)G b (q - p)(l + p 2 )dudqdvdp < C, 



(note that it is important here that h u>q {v,p) = unless u G B 2 ) we get 

(27r/i) n Tr( 7 0(-/i 2 A)0) 
< y x(K,q(v, p)) G b {u - v)G b {q - p)(j)(v + h 2 ab(u - v)) 2 
x (p + h 2 ab(q - p)) 2 dudqdvdp + Cbh 2 . 

From (39) we may now conclude that 

(27r/ i )"Tr(70(-^A)0) < J x{°{ v ,P) +£v(u,q) - C\u\(b~ l + \u\ 2 )) G b (u)G b (q) 

x (p(v + h 2 abu) 2 (p + h 2 abq) 2 dudqdvdp + CWi 2 . (49) 
We now perform the p-integration in (49) and arrive at 

{V(v)+Uu,q)-C\u\(b- 1 + \u\ 2 ))_ ~ 2 
x G b (u)G b (q)<j>(v + h 2 abu) 2 dudqdv + Cbh 2 , (50) 

where we have used that the integral over the term containing q ■ p vanishes and the 
integral over the term containing (h 2 abq) 2 is bounded by h 4 a 2 b < h 2 b. 

We now expand the integrand in (50) in the same way as we did the integrand in 
(35). We finally obtain 

(27r/i) n Tr( 7 0(-/i 2 A)0) < — ^— u n [ \V{v)^ +1 (j){v) 2 dv + Ch^\ 

n + 2 J 

which together with (48) gives (36). □ 
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We shall need the generalization of Theorem 12 and Lemma 14 to a ball of arbitrary 
radius. We also require to know how the error term depends more explicitly on the 
potential. 

Corollary 15 (Rescaled semi-classics). Let n > 3, <fi E Cq +a (R u ) be supported in 
a ball B e of radius £ > 0. Let V G C 3 (B e ) be a real potential. Let H = —h 2 A + V, 
h > and a(u, q) = q 2 + V(u). Then for all h > and f > we have 



Tr[0#0]_ - (27r/i)- n J 0(w)V(w, q)_ dudq 
where the constant C depends only on 



< Ch- n+6/5 f n+4/5 £ n - 6/5 , (51) 



sup \\S a \d a (l)\\ 00 , and sup ||/~ 2 £ |a| 9 a \/|| 00 . (52) 

|a|<ri+4 |q|<3 

Moreover, there exists a density matrix 7 such that 

Tr[0#0 7 ] < {2nh)- n J <P(u) 2 a(u, q)_ dudq + ^-"+6/5^+4/5^6/5 (53) 

and such that its density p y (x) satisfies 

Pj (x) - (2irh)- n LU n \V( X )_\ n/2 < C/,-+9/10 / n-9/10 r 9/10 j (54) 



for (almost) all x G Be and 



J (j){x) 2 Pl {x)dx - (2nh)- n u n J <P(x) 2 \V(x)_\ n/2 dx 



where the constants C > in the above estimates again depend only on the parameters 
in (52). 

Proof. This is a simple rescaling argument. Introducing the unitary operator 
(Uip)(x) = i~ n l 2 i\){l~ l x) we see that <f)H(f) is unitarily equivalent to the operator 

u*<t>H<t>u = f 2 M-h 2 r 2 r 2 A + v^)<i> t , 

where <pe(x) = <j>(£x), and Vf t i(x) = f~ 2 V(£x). Thus 

Tr[0#0]„ = f 2 Tr[M-h 2 r 2 £- 2 A + V u )<j> e ]-. 
Note that <fie and Vf t e are defined in a ball of radius 1 and that for all a 

\\d a Moo = W^d^lU and W^Vf,^ = \\r 2 £ H d a V IU 
It follows from Theorem 12 that 

q)_dudq <Cf 2 {hf~ 1 £- 1 )- n+6/5 , 



Tr[0#0]_ - (2 7 rhf- 1 r 1 )- n J Uv) 2 f<J f , t {u, 



(56) 
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where (Jf,e(u, q) = q 2 — Vfj(u) and where the constant C only depends on the param- 
eters in (52). A simple change of variables gives 

(2irhf- 1 £- 1 )- n J <Pt{u) 2 f 2 a f/ {u,q)_dudq= (2nh)- n J (j)(u) 2 a(u,q)_dudq. 

Thus (51) follows. 

To find the appropriate density matrix 7. We begin with the corresponding density 
matrix 7/^ for (fie(—h 2 f~ 2 £~ 2 A + Vf^fa, i.e. the density matrix, which according to 
Lemma 14 satisfies the three estimates 

Tr[M-h 2 r 2 e- 2 A + V f)t )<l>au] < (27rhf- 1 £- 1 )- n J <p 2 (u)a f/ (u,q).dudq 

+c(hf- 1 r 1 y n+6/5 , 

|p 7/ » - (2nhf- 1 £- 1 )- n u; n \V f , e (x)-\ n/2 \ < W^r 1 )-^ 9 / 10 , 

J cf ) 2 e p lfe -(27rhf- 1 £- 1 )- n u; n J Ux) 2 \V f/ {x)_\ n / 2 dx < C^hf^r 1 )^ 5 . 

The density matrix 7 = U^f^U* whose density is Pj(x) = £~ n p lf e (x/£) then satisfies 
the properties (53-55). □ 



5 Semiclassics for the Thomas- Fermi potential 

We shall consider the semiclassical approximation for a Schrodinger operator with 
the Thomas-Fermi potential V TF (z, r, x), i.e., —h 2 A — V TF . We shall throughout this 
section simply write \^ TF (a;) instead of \^ TF (z, r, a;). Recall that V TF (x) > 0. 

The main result we shall prove here is the Scott correction to the semiclassical 
expansion for this potential. 

Theorem 16 (Scott corrected semiclassics). For allh > and allr 1: . . . , r M e M 3 
with min fc ^ m \r m — r k \ > r > we have 



M 

Tr[-/* 2 A - V TF ]_ - (2^)- 3 J ( P 2 - V^(u))„ dudp - JL g 



(57) 

where C > depends only on zi, . . . , zm, M, and r . Moreover, we can find a density 
matrix 7 such that 

Tr [{-h 2 A - V TF ) 7 ] < Tr [-h 2 A - V TF ] + Ch- 2+1 / 10 , (58) 

and such that 

^ (,/TF)3/2 ) s CTi_5+4/5 (59) 
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and 

J^K^LjJ V^(xf/ 2 dx + Ch-* 1 '*, (60) 
with C depending on the same parameters as before. 

Note that if we choose h = 2~ l l 2 we have from (14) that {Q^h^iy^f 12 = 
p TF /2. The factor 1/2 on the right is due to the fact that we have not included spin 
degeneracy in Theorem 16. 

In order to prove this theorem we shall compare with semiclassics for hydrogen 
like atoms. 

Lemma 17 (Hydrogen comparison). For all h > and all r±, . . . ,r M el 3 with 
mm k ^ m \r m — r k \ > r > we have 

Tr [-h 2 A - V TF (x)] _ - (2vr/i)- 3 J (p 2 - V TF (u))_dudp 

~ Yl ( Tr [- h2A ~ j^Tj + 1 ] - ( 2nh ^ 3 I (p 2 - uT^h + X ) dudp ) 

< Ch- 2+1 / 10 , (61) 
where C > depends only on z±, . . . , zm, M and r . 

The first estimate in Theorem 16 follows from Lemma 17 combined with the exact 
calculations for hydrogen 



Tr 



-tfA - + 1 

\x - r k \ 



-2 



l<n<z k /(2h) 



and 

32tt 2 4 r(7/2)r(l/2) _ z\ 



(2wh)- 3 (p 2 - + l) dudp = 

7 V \u — rh\ / - 



\u-r k \ J- * 15(2tt/i) 3 T(4) 12h 3 
Before giving the proof of Lemma 17 we introduce the function 

M _ L 

t(x) = l(l + J2(\*-r k \ 2 + ®- 1/2 ) (62) 

k=l 

where < £o < 1 is a parameter that we shall choose explicitly in (76) below. Note 
that £ is a smooth function with 

0<£(x)<l, and ||W(aO||oo < 1. 

Note also that in terms of the function d(x) from (18) we have 

|(1 + M)-H Q < |(1 + M(d(x) 2 + tlY 112 )- 1 < e(x) < \{d{x) 2 + tlf 12 . (63) 
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Note in particular that we have 

£{x) > \{l + My 1 m\n{d(x), 1}. (64) 

We fix a localization function <fi G C^(1R 3 ) with support in {\x\ < 1} and such 
that J <p(x) 2 dx = 1. According to Theorem 22 we can find a corresponding family of 
functions 4> u G C(f(R 3 ), u G M 3 , where U is supported in the ball {\x — u\ < £(u)} 
with the properties that 

J <p u (x) 2 £(u)- 3 du = 1 and \\d a (f) u \\ 00 < C£(u)-^, (65) 

for all multi-indices a, where C > depends only on a and 0. 

Moreover, from (21) in Theorem 7 we know that for all u G lR n with rf(-u) > 2£q 
the TF-potential y TF satisfies 

sup \d a V TF (x)\<Cf(u) 2 £(u)~ lal , (66) 

|x-u|<£(u) 

where C > depends only on a, ^i, ... , zm, and M. We have here used the fact that 
if d(u) > 2£ then £{u) < \/hd{u)/A and hence for all x with \x — u\ < £{u) we have 
(note that d(u) < d(x) + \x — u | and v^/4 < 1) 

£{u) < Cd{x) and f(x) < Cf(u). 

Proof of Lemma 17. We note first that we may if necessary assume that h is smaller 
than some constant depending only on the parameters z\, . . . , Zm, M, r . This follows 
from the Lieb-Thirring inequality (7) and the estimate on V TF given in (21) for a = 0. 

In order to control the region far away from all the nuclei we introduce localization 
functions 6_,6 + G C°°(M) such that 

1. e 2 _ + e\ = i, 

2. 0_(i) = 1 if t < 1 and 9_(t) = if t > 2. 
Let 

R = h- 1 ' 2 (67) 

and define $ ± (x) = 6±{d{x)/R). Then <$> 2 _ + $^ = 1. Denote J = (V$_) 2 + (V$+) 2 . 
Then X is supported on a set whose volume is bounded by CR 3 (where as before C 
depends on M) and 

imioo < CR- 2 . 

Using the IMS-formula (10) we find that 
-h 2 A - V TF = $_(-/i 2 A - V TF - h 2 T)<$>_ + $ + (-/i 2 A - V TF - h 2 l)§ + 
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From the Lieb-Thirring inequality the estimates on X and the bound \^ TF (a;) < 
Cd(x)~ A (see (21) with a = 0) we find 

Tr[-/i 2 A - V TF ]- > Tr[$_(-/i 2 A - V TF - /i 2 J)$_]_ - C(/r 3 iT 7 + h 2 R- 2 ). 

On the support of we now use the localization functions <p u . Again using the IMS 
formula (10) we obtain from (65) that 

> J (-h 2 A - V TF - Ch 2 (£(u)~ 2 + R~ 2 )) <p u ^^£{u)~ 3 du. 

We have here used that if the supports of 4> u and </v overlap then \u — u'\ < £(u)+£(u') 
and thus 

£{u) < £{u') + || V£||oo (£(u) + £{u')). 

Therefore, since ||W||oc < 1, we have that £{u) < C£(u') and thus £(u')' 2 < C£(u)' 2 . 
From the variational principle we now get 



Tr[-/i 2 A - V 



TFi 



(68) 



> / Tr[(f) u (-h 2 A-V TF -Ch 2 £(u)- 2 )(f) u }-£(u)- 3 du 

Jd(u)<2R+l 

-C{h~ 3 R' 7 + h 2 R~ 2 ), 

where we have restricted the integral according to the support of $_ and <p u and used 
that, since we may assume that h is so small that R> C, then £(u)~ 2 > CR~ 2 . Note 
that there is no need to write <3>_ on the right, since in general Tr($_A$_)_ > TrA_ 
for any selfadjoint operator A. 

In a very similar manner we get corresponding estimates for the hydrogenic oper- 
ators. In particular, if we choose h so small that R > m&Xk{zk} then on the support 
of $ + we have — z k \x — r fc | _1 + 1 > 0. Thus we have 



Tr 



-h 2 A 



> 



\x - r k 
Tr 



+ 1 



(69) 



d(u)<2_R+l 

Ch 2 R~ 2 . 



-h 2 A- 



Zk 



\x - r k \ 



+ 1 



- Ch 2 £(u)- 2 ^<f) u £{u)~ 3 du 



We shall now get upper bounds similar to (68) and (69). If we again denote by x 
the characteristic function of the interval (— oo, 0] we see from (65) that 

1=1 KX {M-h 2 A - V TF )0 U ) ^ u £{u)- 3 du 

Jd(u)<2R+l 
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defines a density matrix. If we use it as a trial density matrix to get an upper bound 
we obtain 

Tr[-/i 2 A- V TF ]_ < Tr[(-/i 2 A - V TF ) 7 ] 

J Tr[0 M (-/i 2 A-\/ TF )0 u ]_£( M )^. (70) 



d(u)<2R+l 



Similarly 

2a Z k 



Tr 



-h 2 A 



\x - r k \ 



+ 1 



J Tr[0 u (-/i 2 A - — + l)^] J( M )- 3 ^.(71) 



< I Tr 

d(u)<2i?+l 



We now introduce the quantities 

D + {u) := Tr^^-^A-^-C^)- 2 )^^, 



M 

-^Tr[</> U (-/* 2 A 



k=i 



\x - r k \ 



+ 1 



Dju) 



M 

E Tr K~ 



h 2 A- 



Zk 



k=i 



x - r k \ 



+ 1 - Ch 2 £(u) 



and 



-Tr^^-^A-O^]-, 
D sc (u) := (2nhy 3 J <p u (x) 2 (p 2 - V TF (x)),dpdx 



M „ 

-(27r/i)- 3 J] / u (x) 2 (p 
fc=i 



z k 



x - r k \ 



+ 1 



Then from (68), (69), (70), and (71) we have 



M 

Tr [-/ i 2 A-\/ TF ]_-^Tr[-/ i 2 A--^ + lj > / 



D + {u)£{u)- 3 du 

!d(u)<2R+l 

-C(h 2 R- 2 + h- 3 R- 7 ) (72) 



and 

M 



5>[- 



-/i A 



fc=l 



F - r k\ 



+ 1 



Tr[-/i 2 A - T/ TF ]_ > 



D_{u)£{u)~ 3 du 

d(u)<2R+l 

-Ch 2 R~ 2 , (73) 



and from (65) 

/M „ 
(p 2 - \/ TF (x))_dpcfe - (2tt/i)- 3 / (V 



Zk 



\x - r k \ 



1 ) dpdx 



= J D sc (u)£(u)- 3 du. (74) 
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< Ch~ 3 R 



-3 d-7 



(75) 



The same estimates which led to (68) and (69) give 

/ D sc (u)£(u)~ 3 du - I D sc (u)£( U y 3 du 

J Jd(u)<2R+l 

We shall prove the lemma by establishing lower bounds on D + {u) — Dsc(u) and 
D-{u) +D sc (u). 

We consider first u with d(u) < 2£ , where £q is the parameter that occurs in the 
definition (62) of £. We choose 

4 = h, (76) 

where we assume that h is small enough to ensure that £o < 1- In fact, we may also 
assume that £ < r /8. If k is such that d(u) — \u — r k \ we get from (63) that 



\u 



r k \ + £(u) = d{u) + £{u) < d{u) + \{d(u) 2 + £%fl 2 < 4£ < r /2. 



Thus for all x with \x — u\ < £(u), i.e., for all x in the support of <p u we must have 
that \x — r k \ < r /2 and hence \x — Tj\ > r /2. Thus d(x) = \x — r k \ and of the nuclear 
positions r 1; . . . , r M only r k may be contained in the support of <j) u . Since the function 
W k (x) = V TF (x) — z k \x — r k \~ l satisfies the estimate (22) on the support of 4> u we 
have for < e < 1/2 that 



Tr 



<f> u (-(l-e)h 2 A-(l-e) 

z k 



> Tr 

+ Tr <f) u ^ -eh 2 A - e- 

2 a z k 



Z k 



+ (1 



-e)} 



> Tr 



K[ ~h 2 A- 



\x - r k \ 



x-r k 
+ 1 



x - r k 

- (1 - e) - W k {x) - Ch 2 £{u)~ 2 



<Pu 



-Ch- 3 [e£{ufl 2 + e- 3 ' 2 £{uY{l + r 5/2 ) + hh- 3 l 2 £{uY\ 

where in the last line we have used the Lieb-Thirring inequality (7) and the fact that 
Tr[. . .]_ < 0. We therefore have that 



D + (u) > -Ch- 3 [e£(u) 1 ' 2 + e- 3 ' 2 £(uf + hh- 3 ' 2 £(u)- 2 }. 
The quantity D_ (u) is estimated in essentially the same way. We get 



(77) 



Tr 



-h 2 A- 

> Tr 
+Tr 

> Tr 



Zk 



+ 1 - Ch 2 £(u) 



(p., 



\x - r k \ 

l-e)h 2 A-(l-e)V TF )<f> u 
-IrA : ' 



+ (1 - e)W k - Ch< 

F - i k \ 
h 2 A-V^)<p u 

Ch- 3 [e£{ufl 2 + e- 3 ' 2 £{uf + h 5 s- 3 / 2 £(u)- 2 ], 



u 
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and again by the Lieb-Thirring inequality 

TrU u (-/i 2 A - — g + 1 - Ch 2 £(u)- 2 )(j) u ] > -Ch~ 3 [£(u) 3 + hH{u)' 2 } 



Therefore D_(u) satisfies an estimate similar to (77). Since d(u) < 2£ and < 1 we 
have that C _1 £ < £(u) < C-^o- Hence we can replace £{u) by £ i n the above estimates 
if we change the constant C. If we now choose e = i^h 2 (with the choice (76) we 
may assume that h is so small that, indeed, e < 1/2) we get 



D + (u),D_(u) > -ch- 3 [h 2 e~ 1/2 + h-'% /z }. 



-3/»9/2l 



(78) 



By an identical argument using + > a;_ + y_ we get that for -u with rf(-u) < 2£ 
we have 

|Ajc(«)I < Ch- 3 [h% 1/2 + h~HT\. (79) 

We next consider u such that 2£q < <i(w) < 2i2 + 1. We choose again < £ < 1/2 
(not necessarily the same as before) and use the Lieb-Thirring inequality to arrive at 



D+{u) > Tr[0 u (-(l- £ )/i 2 A-y TF )0 u ]_ 

M 

-^Trk(-/* 2 A--^ + l)<d -Ce^ 2 h 



fc=i 



50) 



Since V TF satisfies the estimate (66) on the ball {x \ \x — u\ < £{u)} and 4> u satisfies 
(65) we may use Corollary 15 to conclude that 

Tr[0 u (-(1 - e)h 2 A - V TF ) - (2tt(1 - e) 1/2 /i)" 3 J 2 (x)(p 2 - V TF (x))_dxdp 

< Ch- 3 h 6 / 5 f(u) 19 / 5 £(u) 9 / 5 , 
where C depends only on and M. Hence again using (66) we obtain 

Tr[04-(l- £ )/i 2 A-y TF )0 u ]_ > (2nh)- 3 J ^(p 2 - (x)).dxdp (81) 

-C7r 3 (e/(«) 5 £( M ) 3 + h 6/5 f(u) 19/5 £(u) 9/5 ). 
If d(u) > max^zj.} + 1 then we have that for all k — 1, . . . , M that 

/ <f> u (x) 2 ( ,r - -r-^—r + l) dpdx = 0. 
J V \x-r k \ J- 



Tr 



-/i 2 A- 



+ 1 



and 



On the other hand, if 2£ < d(u) < max^-f^} + 1 then for all k = 1, . . . , M the 
potential zy\x — r^ 1 — 1 satisfies a bound similar to (66) and we may again use 
Corollary 15 to conclude that for all u with d{u) > 2£ we have 



Tr 



-h 2 A 



Zk 



\x - r k \ 



+ 10 



-(2^)- 3 |0 2 (x)(. 



+ 1 ) dxdp 



p - 

\ x ~ r k\ 

< Ch- 3 h & ' b f{u) 19 'H{u) 9 ' b . (82) 
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Hence from (80), (81), and (82) we have for all u with 2£ < d{u) <2R+1 that 

D+(u) - D sc (u) > -Ch^(e- 3/2 h 5 £(u)' 2 + ef(uf£(uf + h 6/5 f(u) l9/ H(u) 9 '' y ) 

= -Ch- 3 (h 2 f(u) 3 £(u) + h 6 / 5 f(u) 19 / 5 £(u) 9 / 5 ), (83) 

where we have chosen e = ch 2 £{u)~ 2 f {u)~ 2 . Note that from the property (64) of £ and 
the definition (19) of / we have 

e < cCh 2 max^M)- 1 , d(uf} < cCh 2 max{4 \ {2R + l) 4 }. 

We see that with the choice of R in (67) and of £ in (76) we may assume that h and 
c are chosen small enough so that e < 1/2. 

In a completely similar way we get for all u with 2£ Q < d(u) < 2R + 1 that 

D_(u) + D sc (u) > -Ch- 3 (h 2 f(u) 3 £(u) + h 6/5 f(u) 19/5 £(u) 9 ^). (84) 

If we now combine (78), (79), (83), and (84) we obtain 

[D±(u) =F D sc {u)]£{u)~ 3 du 



I 



d(u)<2R+l 

-Z\ h 2 ff -l/2 , , -3,9/2 



> -ch-^hX^ + h-x* 

-Ch- 3 J (h 2 f(u) 3 £(u)- 2 + h^ 5 f(u) 19 ^£(u)- 6 ^)du, (85) 



2e <d(u)<2R+l 



where we have used that the volume of the set of u for which d(u) < 2£ is bounded 
by C£\ and that from (63), £(u) > C' 1 ^. Using again the property (64) of £ and the 
definition (19) of / we see that the last integral in (85) is bounded by 

Ch' 3 J h 2 mm{d(u)~ 7/2 , d(u)' 6 } + h 6/5 min{d(uy 31/w , d(u)- 38/5 }du 

2£ () <d(u)<2R+l 

< Ch' 3 J fe 2 min{|u|- 7 / 2 , |w|- 6 } + ^ 6 / 5 min{|w|- 31 /i°, jw]- 38 / 5 }^ 



2£ <\u\ 

< ch-% 1/2 + ch- 9 /% 1/w . 

Thus if we combine (72), (73), (74) ,(75), and (85) we see that the left side of the 
main inequality (61) is bounded by 

C(h~X 1/2 + h- 9 'H~ 1 ' 10 + h-H 9 J 2 + h 2 RT 2 + Ch- 3 R- T ) 

< C {h- 3 ' 2 + h- 19 / 10 + h~ 3 ' 2 + h 3 + h 1 ' 2 ) < ch- 2+1 / 10 . 

Note that the choice of £ has not been optimized. □ 
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Proof of Theorem 16. As mentioned just after the statement of Lemma 17 the first 
estimate in the theorem is a consequence of the lemma. It remains to prove the 
existence of a density matrix 7 with the stated properties. 

We note first that we may as before, if necessary, assume that h is smaller than 
some constant depending only on the parameters Zi, . . . ,zm, M, and r . Otherwise 
we simply choose 7 = 0. That this is an acceptable choice follows from the Lieb- 
Thirring inequality (7) and the fact that the estimate (21) for a = implies that 

d(V tf ) 3/2 ) <C. 

To construct 7 we shall again use the localization family (j) u with the properties 
given in (65). As in the previous lemma we shall choose R = hr 1 ! 2 and £ = h 
(although it is not a requirement that they should be as before) . We shall choose 7 of 
the form 

7 = / <f) u i u <j) u £{u)~ 3 du, (86) 

Jd{u)<R 

where 7„ is a family of density matrices which we shall now choose. Note that the 
first condition in (65) implies that 7 is then a density matrix. 

If 2£ < d{u) it follows from (65), (66), and Corollary 15 that we may choose 7„ 
such that (53), (54), and (55) hold when V = V TF , = <j> u , £ = £{u) and / = f(u). 

If d(u) < 2£ we simply choose 

In = X (<Pu (-h 2 A - V TF ) <t> u ) , 

where x is again the characteristic function of the interval (— 00, 0]. I.e., 7„ is the 
projection onto the non-positive spectrum of 4> u (— h 2 A — V TF ) <f> u . Here we are con- 
sidering; oL ClS db multiplication operator. 

We shall first prove that for d(u) < 2£ we have 

/ (4fa.f S <Ch-%> 2 . (87) 

From the Lieb-Thirring inequality (7) and the estimate (21) with a = 0, we conclude 
that 

> Tr (0 U (-h 2 A - V TF ) <p ulu ) > |Tr (-/i 2 A(0 u7u M )) - ChrHf \ 

where we have used that d(u) < 2£ implies that £(u) < C£ - The density of the 
operator 0„7 U U is (fi 2 p-y u ■ Thus, using the Lieb-Thirring inequality in the formulation 
(6) we arrive at (87). 

Using (87), Holder's inequality, the support property for <p u , (63), and the second 
property in (65) we arrive at the estimates 

||0^ 7 Jli<Crt(«) 3/2 and ||0^ 7 J| 6 /5 < ChrH(u\ (88) 
for u with d(u) < 2£ . For these u we also have 

H0^ TF ) 3/2 ||i<^f and ||0^ TF ) 3 / 2 || 6/5 <^ O , (89) 
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where we have used that from (21) with a = 0, V TF (x) < Cd{x)~ 1 . 

We are now ready to prove that the density matrix 7 has the stated properties. 
We begin with proving (60). The density of 7 is 



<fru(x) 2 p lu (x)£(u) 3 du. 



d(u)<R 



From (55) we see that for d(u) > 2£ we have 



(xfp lu (x)dx < 



67T 2 h 3 

and from (88) and (89) we get for d(u) < 2£ that 

1 



u (x) 2 V TF (x) 3 / 2 dx + Ch- 2+1 ' b f{ufH{uf^ 



^(x) p^ u (x)dx < 



Q7T 2 h 3 



u {xfV T¥ {xfl 2 dx + Ch-Hl' 2 . 



Hence using the first property of (f> u in (65) we obtain 
J Pl (x)dx < J V TF (xf 2 dx + ChT 3 

+C 
1 



-0 



(u) 3 du 



d(u)<2£ 

h- 2+1 ^f(uf 5 £(u)~^du 



< 



6n 2 h 3 

+ch- 2+1 / 5 



2£ <d(u)<R 

v TF (x) 3 / 2 dx + ch~ 3 e 3 /2 



mm{d( U y 21/10 ,d( U y 18/5 }du 



2£ <d(u)<R 



where we have inserted the choice £ = h, used that h is small, and controlled the 
integral over the region 2£ < d(u) < R in a way similar to the integral in (85), using 
the properties (19) and (64) of / and £. 

We now come to the proof of (59). If we use the Hardy-Littlewood-Sobolev in- 
equality (9) we see that it is enough to estimate the 6/5 norm 



Pi 



Qn 2 h 3 



TF\3/2 



< 



6/5 



d(u)<R 



K Plu - 



1 



Qn 2 h 3 



TF\3/2 



£{u)- 3 du 



+ 



J 

JdO 



Qn 2 h 3 



6/5 

u)~ 3 du. 



6/5 



>d{u)>R 

If we use (88) and (89) when u with d(u) < 2£ , (54) when 2£ < d(u) < R, and (66) 
when d(u) > R we obtain 

1 



Pi 



6n 2 h 3 



TF\3/2 



6/5 



< Ch' 3 £ + Ch' 21 / 10 J f{u) 2l / w £{u)~ 7 ' b d 

2£ <d{u)<R 

+Ch' 3 [ f(u) 3 £(uy l ' 2 du. 

Jd(u)>R 
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Using as before the properties (19) and (64) we see that the first integral above is 
bounded by a constant and the second integral (since we may assume that R > 1) is 
bounded by R~ 3 . Thus using the Hardy-Littlewood-Sobolev inequality (9) we obtain 



D (* - ^' TF > 3/2 ) £ c 



1 



P-y 



TF\3/2 



<C(h- 5+4 / 5 + h- 6 R- 6 ). 



6/5 



Finally we turn to proving (58). From the definition of 7, (51) and (53) we obtain 
Tr [{-h 2 A - U TF ) 7 ] = / Tr[0 M7u tt (-^A-\/ TF )]£H- 3 ^ 

< f Tr [4> u {-h 2 A - V TF ) 0j _ £(u)- 3 du 

Jd(u)<R 

+ Ch -2+i/5 f f( u ) 19/5 £(u)- 6/5 du. (90) 

J2t <d(u)<R 

On the other hand, from (68) (used with 2R + 1 replaced by R) we get 

Tr[-/i 2 A - U TF ]_ > / Tr[0 M (-/i 2 A- V rF - Ch 2 £(u)' 2 ) ^^(u^du 

Jd(u)<R 

-C(h 2 R- 2 + h- 3 R- 7 ). (91) 
Appealing to the Lieb-Thirring inequality (7) we see that for all < e < 1 we have 

Trl^ (-h 2 A - V TF -ChHiu)- 2 )^ > (l-e)Tr[<f> u (-h 2 A-V TF )<f> u }_ 

- Ce-^ 2 h 2 £(u)- 2 

-Ceh' 3 [ (V TF (x)f 2 d. 
J\x-u\<e(u) 



X. 



Using (21) and (66) both with a = we find that 



V TF (xf 2 dx < 



CI 



1/2 
5 



if d(u) < 2£ 



<\ x -u\<e(u) l Cf(u) 5 £(uf, ifd(u)>2£ ' 

If d(u) < 2£ we choose e = d^h 2 (as we did just before (78)) and we get 
Ti[<p u {-h 2 A-V^ F -Ch 2 £{u)' 2 )<p u }_ > Tr[0 u (-h 2 A - U TF ) 0, 



Cbr v i 



1,-1/2 



(92) 



If d{u) > 2£ we choose e = h £{u) f(u) (as we did just after (83)) and we get 

Tri<p u {-h 2 A-V TF -Ch 2 ^)- 2 )^ > Tr[0 u (-/i 2 A-U TF )0 u ]_ 

- Ch- l f(uf£(u) (93) 
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Thus from (90), (91), (92), and (93) we get 
Tr [(-h 2 A - V TF h] < Ti[-h 2 A-V T¥ ]^ + C{h- l r 1,2 + h 2 R' 2 + h~ z R^) 

+C j h- 2+1 / 5 f(u) 19 / 5 £(u)- 6 / 5 + h- 1 f{u)H{u)- 2 du. 

d(u)<R 

The estimate (58) now follows from a calculation almost identical to the one given 
right after (85). 

□ 

6 Proof of the Scott correction for the molecular 
ground state energy 

The proof of the Scott correction Theorem 1 is now a fairly standard application of the 
results in the previous sections. We begin with giving the proof of the lower bound. 

Lemma 18 (Lower bound). Let R and Z be as in the statement of Theorem 1. 
Then, the ground state energy for a neutral molecule satisfies 

E(Z, R) > E TF (Z, R) + lJ2 Z l + 0{\Z\ 2 - 1 '™). 

j 

Proof. The starting point is the Lieb-Oxford inequality (11), from which we conclude 
that if if) is a Z-particle (N = Z) wave function we have 



H(Z, R)fj > (*!>, [ - - V(Z, R, Xl )}^) + D(p^) -C j pf . 
i=i J 



In order to bound the last term we use the many-body version of the Lieb-Thirring 
inequality (8). For all < e < 1/2 we have 



^ej^-\^)- C \pf > efpf-cfptl 

i=l J J 



4/3 

1/2 / r \ 1/2 



> -e- l C J Pi> = -Ce- x Z. 

Here we have used Holder's inequality for the p 4//3 integral and used that ip is a 
Z-particle state. Thus 

(i>,H(Z,R)i>) > /^,^(-(l-e)lA i -y(Z,i2,x0)^\+£>W)-C'e- 1 Z 
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-D(p TF (Z, R,-))-Ce~ l Z 
> 2 Tr[ - 1 (1 - e)A - V TF (Z, R, •)]_ - D(p TF (Z, i?, •)) - Ce _1 Z. 

Here we have applied (13), the fact that the Coulomb kernel is positive definite such 
that D(p — p TF ) > 0, and the Fermionic property of the wave function. 
If we now use the scaling property (15) we find that 

Tr[ - 1(1 - s)A - V TF (Z, R, •)]_ = \Z\^Tr[ - 1 (1 - e)\Z\-^/\ - V TF (z, r, •)]_, 

where z = (zi,...,zm) and r = (r 1? ...,r M ). Using now (57) (with h = 
v/(l - e)/2\Z\- 1 / 3 ) and (27) we see that 

2Tr[ - 1(1 - e)A - V TF (Z, R, •)]_ = (1 - e)- 3 ' 2 \Z\ 7 / 3 {E TF (z, r) + D(p TF (z, r, •)) 

| 7 |2 M 

+d-s)- il -^-j:4+o(\zr^) 

= (l-er^(E^ F (Z,R) + D(p^ F (Z,R r )) 

M 

+\{l-e)-^ Z l + 0{\Z\ 2 - llm )- 

k=l 

We have here used the TF scaling E TF (Z, R) = \Z\ 7 / 3 E TF (z, r) and D(p TF (Z, R, •) = 
\Z\ 7 ^D (p TF (z, r, ■)). Choosing e = \Z\~ 2 I 3 completes the proof of the lemma. □ 

Lemma 19 (Upper bound). Let R and Z satisfy the conditions from Theorem 1. 
Then, the ground state energy for a neutral molecule satisfies 

E(Z, R) < E TF (Z, R) + \ Z) + OQZl 2 - 1 / 30 ). 

3 

Proof. The starting point now is Lieb's variational principle, Theorem 5. By a simple 
rescaling the variational principle states that for any density matrix 7 on L 2 (R 3 ) with 
2Tr7 < Z we have 

E(Z,R) < \Zf 3 (2Tr [{-\\Z\- 2 ' 3 /\ - V{z,r,x)) 7] + \Z\D(2\Z\~ l Pl )) . 

As for the lower bound we bring the TF-potential into play 

\Z\- 4/3 E(Z,R) < 2Tr [(-||Z|- 2/3 A - V(z,t,x))j] + \Z\D{2\Z\~ 1 p 1 ) 
= 2Tr[(-||Z|- 2 / 3 A-y TF (z,r,x)) 7 ] 

+\Z\D (2|Z|-V 7 - p TF (z, r, •)) - \Z\D(p TF (z, r, •))• (94) 
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We now choose a density matrix 7 according to Theorem 16 with h = \J\j1\Z\ x / 3 . 
Note that with this choice of h we have that 

(67r 2 /r 3 )-V TF (z,r,:r) 3 / 2 = \Z\p TF (z, r, x)/2. 

Since J p TF (z,r,x) = Y^jLi z i = 1 we see fr° m (60) that 

2Tr7 < \Z\ + C|Z| 2/3 - 1/15 = |Z|(1 + C\Z\' 1/3 ~ 1/15 ). 

Thus if we define 7 = (1 + C^- 1 / 3 - 1 / 15 )- 1 7 we see that the condition 2Tr7 < \Z\ is 
satisfied. 

Using (59) we see that 

\Z\D (2|Z|- Vy - p TF (z, r, •)) < C\Z\ 2 ^ 4 ' l \ 

and thus 

\Z\D (2\Z\- 1 p 1 -p TF (z,r,-)) 
< C|Z|(1 + ClZr 1 / 3 " 1 / 15 )- 2 /} (2|Z|-^ - p TF (z, r, •)) 

+C | Z |l/3-2/l 5jD ( p TF (z? r < CIZI 2 / 3 - 4 / 15 , (95) 

where we have used the triangle inequality for a/Z>, and that D (p TF (z, r, •)) < C. 
Finally, if we use (57), (58), and (27) we arrive at 

2Tr[(-i|Z|- 2 / 3 A-\/ TF (z,r,x))7] < \Z\ (E TF (z, r) + JJ(p TF (z, r, •)) 

1712/3 M 

+^-E^ 2 +°(^i 2/3 " /30 )- 
k=i 

Since E TF (z, r) < C and D (p TF (z, r, •)) < C we see that the same estimate holds for 
7 replaced by 7. If we insert this estimate together with (95) into (94) and use again 
that E TF (Z,R) = \Z\ 7 / 3 E TF (z, r) we arrive at the upper bound in the lemma. □ 

A Appendix: Results on the new coherent states 

Before we prove the trace formula (10) and the representation (11) we need some 
simple lemmas. The proof of the first one is a straightforward calculation which it is 
left to the reader to check. 

Lemma 20. Let Q u q be defined as in (28), then its integral kernel is 

g Utq (x,y) = (nh)- n l 2 e-<^- u ) 2+iq{x - y)/h -^ {x - y)2 . (96) 
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Lemma 21. Let B e R, -81,2 G ^ n arad A = _B + -B]£ — ihB 2 V be a linear combina- 
tion of the identity, multiplication and momentum operator. Let f be a polynomially 
bounded, measurable function on M, n with values in ffiL Then, 



f(B + B 1 x-ihB 2 V)(x,y) = J f(B + B 1 (^)+B 2 p)e l 



(x-y)/h_ 



dp 



{2nh) n 
as a distributional kernel. 

Proof. First assume that B 1 and B 2 are not orthogonal. By applying the unitary 
transformation (Uip)(x) = exp [— 2h ^ B ^ (Bix) 2 ]ip(x) and the Spectral Theorem we 
get U- 1 f(B + B±x - ihB 2 V) U = f{Bo~ihB 2 V). Hence, 



f(Bo + B lX -ihB 2 V)(x,y) 

f(B + B 2 p) e^^BT^ 



/ 

J f(B + B 1 ( x -¥)+B 2 p) e^ 



((Bix) 2 -(Bij/) 2 )+ip(x-j/)//i dp 



(2nhy 

ip(x-y)/h d P 



{2Tth) n ' 

The case B\ ■ B 2 = follows then, say by continuity. □ 

Proof of Theorem 9. The proof is actually a fairly standard exercise in calculations 
with Fourier integrals, but we shall do it here carefully. 

Let us for simplicity call Q = G u ,q- Since / is polynomially bounded we have that 
(1 + A 2N )~ l f(A) extends to a bounded operator when iV is a large enough integer. 
From the explicit expression (96) for the integral kernel of Q we immediately see that 
(l + A 2N )g(x,y) is in L 2 (W n x W n ). Thus, the operator {l + A 2N )g is Hilbert-Schmidt. 
It follows that f{A)Q is Hilbert-Schmidt. Likewise (1 + |a;|) -JW V(:r) is bounded for M 
large enough, and thus QV(x) extends to a Hilbert-Schmidt operator. Moreover if we 
define 

f e (s) = f\s)e- £s2 and V £ {x) = V(x)e- £X \ 



we have 



limTr 



Gf £ (A)QV £ \ = Tr [Gf{A)QV . 
The trace on the left can be immediately calculated from Lemmas 20 and 21. □ 



Proof of Theorem 10. We proceed as in the previous proof. We consider first the 
case V = 0. We have that (1 — h 2 A)~ 1 F(—ih'V) is a bounded operator and thus 
(f>F(—ihV)(pG is a Hilbert-Schmidt operator. It follows moreover that if we define 

f £ (s) = f(s)e- £s2 and F £ (p) = F{p)e~^\ 



then the operators 



f £ (A)G and 4>F £ {-ihV)4>G 
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converge in Hilbert-Schmidt norm as e — > 0. We therefore have 

Tr[£/(A)£0FHW)0] = ]im]imTr[gfs(A)g<j>F e (-ihV)</>]. 

The trace on the right can be written as an absolutely convergent integral 

(2nh) 2n Ti[g f s (A) Q <f>F £ (-ihV)<f>] 

= I fs(B + B 1 ^ + B 2 p)e ip{ - x ~ v ^ h+i ^ z ~ z '^ h g Uyq {y, z)g uA {z'x) 

x <p(z)(p(z')F £ (ri) drjdxdydzdz'dp 

We now perform the integration in the variable x — y and we rename the integration 
variable (x+y)/2 as v. We subsequently change variables so that (z, z', rf) are replaced 
by (z + h 2 ab(u — v) + v, z' + h 2 ab(u — v ) + v , rj + p + h 2 (ab)(q — p)). This eventually 
gives 

(27r/i) 2 >/4&)"/ 2 Tr[£ fs(A) Q 0F £ (-iW)0] 
= J fs(B + B lV + B 2 p)F £ (r] + w(p, q )) e - b ^ 2 - b ^ 2 (97) 

x 4>{z + w(v, u))(f)(z' + w(v, u))e iri{z - z ' )/h e~ b ^) -jk {z+z '^ dzdz'di]dvdp, 

where to avoid lengthy expressions we have introduced the function w(s,t) = s + 
h 2 ab{t — s) for s,t G M. n . We now expand F(w(p, q) + rf) around w(p, q). I.e., we write 



F £ (w(p, q) + r})= F(w(p, q)) + r}- VF(w(p, q)) + R F (w(p, q),r}) 



e -s(w(p,q)+Tj) 2 



where according to the assumption that all second derivatives of F are bounded we 
have that the remainder term satisfies \Rp(p + h 2 ab(q — p) , r])\ < Ci] 2 , with a constant 
depending on the bound on the second derivatives of F. We shall estimate the error 
coming from the remainder term below, but we first consider the contribution from 
the two main terms. If we insert the two main terms in the expansion into the integral 
(97) above we see that they give integrals in which the 77 integration can be performed 
explicitly. After letting (5,e^0we obtain for these leading terms 

{2ixh) n J f(B + B lV + B 2 p)F(w(p, q )) e -b(p- q ) 2 -b(u-vr 

4>(z + w(v, u)) 2 e~~h?t> z2 dzdvdp. 

Here we may now expand the function <fi 2 

(f)(z + w(v, u)) 2 = (j)(w(v, u)) 2 + z ■ V<f)(w(v, u)) 2 + R<f,{w(v, u),z), 

where \R^{h 2 ab{u — v) + v, z)\ < Cz 2 with a constant that depends on the bound on 
0, and its first and second derivatives. Hence the integral in (98) may be written as 

(27r/i)"(6/47r) n / 2 J f(B + B x v + s 2 p) e -6(P-9) 2 -K«-^) 2 F(w(p, q)) 
x [<f)(w(v , u)) 2 + Ei(u, v )] dvdp, 
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where \E 1 (u,v)\ < Ch 2 b. 

We now return to estimating the contribution from the remainder term R F . Note 
that for all integers 1 < k < (n/2) + 2 we have 



(l + V 2 /(bh 2 )) k 



J <f>(z + w(v,u))i>(z' + w{v ) u))e ir i {z - z '^ h e~ b ^) d(z - z') 



< C J (l-b- 1 A g . g ,) k (^{z + w(v,u))cj)(z' + w(v,u))e 6 M ^ 



d(z - z') 



< Cb- n/2 , 



where we have used that b > 1. Here the constant depends on the bound on 0, and 
its first n + 4 derivatives. Thus the relevant contribution to the integral (97) coming 
from the error term Rp can be estimated by 

R F (p + h 2 ab(q - p),r])(f)(z + w(v, u))(f>{z' + w(v, u))e l ^ z - z ' ),h 

x e~ b ^) ~^ {z+z ' )2 dzdz'dr] 

< Cb~ n/2 J(l + r] 2 /(bh 2 ))~ k 7] 2 e-^ 2 dzdr] < Cb n/2 h 2n bh 2 , 

where we have chosen k so as to make the integral finite. We can always do this 
without violating 1 < k < (n/2) + 2. Thus after taking the limit 5 — > we obtain the 
statement of the theorem. 

The case when F = and V ^ is similar but much simpler since we may start 
with Theorem 9 and expand V. □ 

Proof of theorem 11. Since a is a sum of a function of q and a function of u it is 
enough to consider only one of the terms, say, V. Let as before G b (x) = (b/-n) n / 2 e" bx . 
It follows immediately from (96) that 



2 dq 



J Gu,q(& - u)G UtQ 

As a consequence we have 



G b (x - u), 

)" 

d( l M L 2 



(2tt/i)' 



= (1 — h ab)(x — u)G b (x — u). 



(99) 



V(x)- J g u , q (v(u) + l b AV(u) + VV(u)-(x-u))g u , q -^ 

= jG b (x-u)(y(x)-(y(u) + l b AV(u) 
+ (1 - h 2 ab)VV(u) ■ (x - u)^du. 
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Using Taylors' formula we have 

V(x) =V(u) + W(u) -(x-u) + U2 9idjV(x)(xi - u t ) (x, -Uj)- K x (x, u) 



where 



Ki(x,u) = \ I ^didjd k V(u + t(x - u))(xi -Ui)(xj -Uj)(x k -u k )(l - (1 -t) 2 )dt. 
Jo ij,k 

Since f XiXjGb(x)dx = we have 

V(x) - J G u , q {V{u) + ±AV{u) + W(u) ■ (x - «)) Gu, q -^ 
= J G b (x -u)(±(AV(x) - AV{u)) + h 2 abW{u) ■ (x - u) - TZ^x, ufjdu 
= J G b {x - u) (i (AV(x) - AV(u)) - \h 2 aAV{u) - K^x, «)) du 

where the last identity follows by integration by parts. The theorem now follows easily 
since a < b and 



AV(x) - AV(u) 



IX 



diAV(u + t(x — u))(xi — Ui)dt. 



□ 



B Appendix: A localization theorem 

Theorem 22. Consider <p G C£°(IR n ) with support in the ball {\x\ < 1} and satisfying 
f 4> 2 (x)dx = 1. Assume that £ : R™ — > R is a C 1 map satisfying < £(u) < 1 and 
||V£||oo < 1- Let J(x,u) be the Jacobian of the map u \— > jr^-, i.e. 



det 



Xj - Ui)dj£(u) 
£(u) 



J(x,u) = £(u)- n 

We set (f) u (x) := <^(fgf) y/ J(x,u)i(u) n ' 2 . Then, for allxeM 

[ (j) 2 u {x)£{u)- n du = 1 

and for all multi-indices a we have 

W^uWoo < £(u)~\ a \C a max II^IU, 

where C a depends only on a. 



(100) 



(101) 
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Proof. In order to prove (100) it is of course enough to consider the case x — 0. 
The identity follows from the change of variables formula if we can show that the 
map F : R™ — > R n given by F(u) = —u/£(u) is a bijection of F^ 1 ({\x\ < 1}) onto 
{\x\ < 1}. 

The map is, in fact, onto R n since F(0) = and > \u\. Hence for all «gK™ 

there exists t G R with — 1 < t < such that F(tu) = u. 

That the map is also injective on F~ l {{\x\ < 1}) follows since for u ^ we may 
write F(tu) = —g(t)u and the map g : R — > K. is monotone increasing for all t for 
which |(yf(t)||-u| < I. In fact, g(t) = t/£(tu) and thus 

g'(t) = l{tu)~ l - U(tu)- 2 V£(tu) ■ u = l{tu)- 1 [1 - tl{tu)- 1 V£(tu) ■ u] 
> e(tu)- 1 [l-\\V£\\ oo \g(t)\\u\]>0. 

Note that <f> u (x) = 4> u (f^y), where 

4> u (x) = (f)(x) det [xidj£(u) + 5^ . 
The estimates (101) follow since Hc^uHoo < C max.\p\<\ a \ H^^Hoq. □ 
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